Diffusive transport and self-consistent dynamics in coupled maps 
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The study of diffusion in Hamiltonian systems has been 
a problem of interest for a number of years. In this paper 
we explore the influence of self-consistency on the diffusion 
properties of systems described by coupled symplectic maps. 
Self-consistency, i.e. the back-influence of the transported 
quantity on the velocity field of the driving flow, despite of 
its critical importance, is usually overlooked in the descrip- 
tion of realistic systems, for example in plasma physics. We 
propose a class of self-consistent models consisting of an en- 
semble of maps globally coupled through a mean field. De- 
pending on the kind of coupling, two different general types of 
self-consistent maps are considered: maps coupled to the field 
only through the phase, and fully coupled maps, i.e. through 
the phase and the amplitude of the external field. The analo- 
gies and differences of the diffusion properties of these two 
kinds of maps are discussed in detail. 
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I. INTRODUCTION 

Understanding transport is a problem of considerable 
practical and theoretical interest in a great variety of 
fields ranging from geophysics to chemical engineering 
and plasma physics. In some situations one can safely 
consider the simple case of passive transport in which 
the transported quantity does not affect the advecting 
flow [1,2]. In the case of a scalar passive field 0(x, t), the 
evolution equation is the advection-diffusion equation 
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(v0) = A)V 2 6, 



(i) 



where v(x, t) is the velocity field, Dq the molecular dif- 
fusivity, and O represents the scalar concentration, e.g. 
the temperature of the fluid or the concentration of a 
pollutant. The domain of applicability of Eq. (1) is lim- 
ited by two important physical assumptions: has to 
be inert (no possible chemical or biological reactions are 
considered) and passive (there is no feedback on the ve- 
locity). Reactive processes can be taken into account 
by adding to the right hand side of Eq. (1) a function 



/(0) modeling the reaction kinetics [3]. This leads to 
the so-called advection-reaction-diffusion (ARD) equa- 
tions widely used in the modeling of chemical and bi- 
ological systems including combustion, diluted chemical 
reactions, and population dynamics [4]. 

Taking into account the feedback of on v, i.e., the 
problem of active transport, is in general more compli- 
cated as this involves the equation of motion for v. Be- 
cause of this, active transport is also called self-consistent 
transport. A well-known example is fluid convection in 
the Boussinesq approximation [5]. In this case, repre- 
sents the fluid temperature which is an active scalar in 
the sense that it modifies the velocity field through the 
buoyancy force in the Navier-Stokes equation for v. 

The goal of the present paper is to study the problem of 
diffusion of active scalars. In particular, we are interested 
in the relationship between self-consistent chaos and dif- 
fusion due to chaotic advection. The study of diffusion re- 
quires an accurate numerical integration of the equations 
of motion for very long times and many initial conditions. 
A common strategy to bypass this technical difficulty is 
to describe the time-continuous equations of motion with 
a discrete-time map. Here we follow this approach and 
study the problem of diffusion in self-consistent symplec- 
tic maps. In the remaining of this introduction we discuss 
in some detail the problem of self-consistent transport in 
fluids and plasmas. The intent of this discussion is to 
provide a physical motivation for the use of globally cou- 
pled maps for studying self-consistent transport. 

One of the simplest physical examples of active trans- 
port is two-dimensional incompressible flows. This mo- 
tion is described by Navier-Stokes equation (1) in which 
represents the vorticity, ( = V x v. Plasma physics 
is another area in which the problem of self-consistent 
transport is crucial. For example, in the Vlasov descrip- 
tion of an electron plasma [6] (in a uniform neutralizing 
ion background) the system is described by the phase 
space electron distribution function / which, for a one- 
dimensional system, evolves according to the Liouville 
equation 



d t f + ud x f + d x <t>d u f = D dlf, 



(2) 
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where the term on the right hand side is a Fokker-Plank 
collision operator, and (a;, u) are the phase space coor- 
dinates. This equation is analogous to Eq. (1) if one 
identifies with /, and v with the transport velocity 
in phase space (u,d x <j>). In this case the self-consistent 
coupling is provided by the Poisson equation 

a l<t> = J f{x,u,t)du - 1 , (3) 

where the right hand side is the charge distribution in- 
cluding the fixed neutralizing ion background. That is, 
the dynamics of an electron plasma is an active transport 
problem in which the transport velocity of the distribu- 
tion function / in phase space is determined by / through 
Poisson's equation. 

The previous ideas on self-consistent transport can be 
reformulated within the Lagrangian description accord- 
ing to which transport is described in terms of individual 
particle trajectories instead of scalar fields and distribu- 
tion functions. The Lagrangian description is important 
because it is the natural description to formulate the self- 
consistent transport problem in terms of symplectic maps 
which are the main objects of study in the present paper. 
As it is well-known [7], the Lagrangian formulation of (1) 
is the Langevin equation 

^=v(x,i) + v /2A^), (4) 

describing the motion of a test particle (the tracer), 
where rj is a normalized Gaussian white noise with zero 
mean and delta correlated in time: 



shown below, the mean field will depend on the mean 
distribution of particles only, thus (6) formally reduces 
to 

—j^ = v ext (x„) + v(x„ - X), 

(7) 

^=F(X,{x,}), 

where we have included the possible contribution of an 
external field v eKt . 

Recently, a mean-field description of this sort has been 
proposed to study self-consistent transport in fluids and 
plasmas [10-12]. In the plasma physics context, this 
approximation known as the single- wave-model (SWM), 
consists of simplifying the self-consistent coupling be- 
tween / and 4> given by Poisson's equation (3). 

The SWM is a general model for the description of 
marginally stable fluids and plasmas [13,14]. Also the 
model bears many interesting analogies with coupled os- 
cillator models used in statistical mechanics [15]. As 
such, it is an insightful model to explore the problem 
of self-consistent chaos, and will be our starting point 
for the construction of the self-consistent symplectic map 
models in the present paper. 

The transition to the symplectic map description is 
eased by first writing the SWM as a full TV + 1 Hamilto- 
nian system in the particle coordinates (xj,Uj) and the 
mean field degrees of freedom [16,10] 

dx k _ cTH_ dpk _ ,o\ 

dt dpk ' dt dxk ' 



< Vi(t)Vj(t') >=5 tJ 5(t-t r ). 



(5) 



The passive scalar nature of 9 in Eq. (1) reflects in the 
absence of coupling in the Lagrangian equations of mo- 
tion (4). However, in an active transport problem, the 
self-consistent coupling between the field and the trans- 
port velocity leads to a nonlinear coupling between the 
Lagrangian equations of motion. The fact that particles 
interact (and usually through long-range interaction) im- 
plies that the phase space evolution of particle n, depends 
on the position of all the TV particles 



dx„(t) 
dt 



= v(xi(i), ...,x N (t)), 



(6) 



and therefore the system has a phase space of dimen- 
sion proportional to TV. This is the well-known TV-body 
problem that arises in many fields of physics, including 
gravitational systems in Astronomy [8] , point vortices in 
two-dimensional fluid dynamics [9], and atomic physics. 

An approximation of the TV-body problem (6), which 
is often used, is a mean field type approximation in which 
the interaction among particles occurs through a global 
variable X function of all the particles. In the examples 



dO _ dU 
dt~~dJ' 
with Hamiltonian 



dJ_ 
~dt 



dU 
~d6 



N 



i=i 



1 
2T~ 



pj -2T 3 



COs(Xj 



(9) 



nj . (io) 



From (10) it is clear that the SWM model consists of TV 
pendulum Hamiltonians mean-field-coupled through the 
amplitude J and the phase 6. Therefore, being the stan- 
dard map the symplectic discretization of the pendulum 
Hamiltonian, the models studied here will consist of an 
ensemble of TV standard maps. In the absence of cou- 
pling, i.e. ignoring self-consistency, J and 6 would be 
constant and the parameters of the standard map would 
be fixed numbers. However, when self-consistency is in- 
corporated, J and 6 become dynamical variables (also 
described by symplectic maps) , and this leads to a depen- 
dence of the parameters of the maps on {xi,x 2 , . . .x N } 
which gives rise to a global coupling of the maps. The 
specific form of this coupling will be discussed in Sec- 
tion III, where we present a systematic discussion of the 
map models in terms of generating functions. 
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The remaining of this paper is organized as follows. 
In Section II we briefly review the diffusion properties in 
the case of the passive scalar, with particular emphasis 
on the standard map. As mentioned before, in Section 
III we introduce the two self-consistent systems studied 
in this work. Sections IV and V are devoted to the dis- 
cussion of the numerical results. Section VI contains the 
conclusions. 



II. A BRIEF REVIEW OF THE DIFFUSION 
PROPERTIES OF PASSIVE SCALARS 

There exists a huge literature about the transport 
properties in the passive scalar limit [1,2,17]. On the 
contrary, there are very few attempts in the study of 
the self-consistent diffusion. The aim of this section is 
to recall the main results on the diffusion problem for 
passive scalars in order to compare them with the self- 
consistent diffusion that will be considered in the next 
sections. It is remarkable that the Lagrangian motion 
can exhibit nontrivial behavior even for a very simple 
velocity field v(x, t) [17,18]. Complex behavior can be 
originated both from chaotic advection (which is in gen- 
eral possible for stationary 3d flows or for time-dependent 
2d flows) and/or from combined effects of the molecular 
diffusivity and the advection velocity. Under very gen- 
eral conditions (see below) the large scale field (9), which 
is obtained as an average of a field G evolving through 
Eq. (1) on a volume whose dimensions are much larger 
than the typical length scale of v, obeys at large times a 
diffusion equation 



dt 



d 2 (9) 
dxidxj 



(11) 



with eddy diffusivity Dfj. In other words, the effect of 
the velocity field at large scales and time is the renor- 
malization of the transport coefficient Dq. It is easy to 
understand the origin of (11) in the Lagrangian frame- 
work. Starting from (4), taking the average over many 
tracers, one has 



which defines the Lagrangian correlation time Tl- From 
(12) one recovers the Taylor result [19] 

(( Xi (t) - x,(0)) 2 ) « 2 (D + (vf)T L ) t = 2D?t, (14) 

which defines the eddy diffusivity in Eq. (11). 

Beyond the above typical scenario one can have anoma- 
lous dispersion, i.e. 



{( Xi (t) - Xi (0)Y 



t 



2v 



(15) 



with v 7^ 1/2. The case v < 1/2, called subdiffusion. 
Superdiffusion {y > 1/2) has been observed in incom- 
pressible flows [21], random shear flows and, as we will 
see, also symplectic maps [20,22]. Anomalous diffusion 
can occur only if some of the hypothesis of the above ar- 
gument breaks down. Practically, this can be due to two 
different mechanisms: 

(a) Infinite variance of the velocity: (v 2 ) = oo. 

(b) Lack of decorrelation: Tl = oo. 

The first condition, which leads to the class of Levy 
flights, is not particularly realistic in physical systems, 
because it requires infinite energy. We will not discuss 
here this behavior. Case (b) is physically more relevant. 
From (13) one sees that anomalous superdiffusion is pos- 
sible only if Cu(t) goes to zero slower than t -1 . Un- 
fortunately, the behavior of CV, (t) is generated by the 
Lagrangian dynamics itself so it is not trivial at all, in 
the absence of molecular diffusivity, to determine whether 
the diffusion process will be standard or anomalous. If 
molecular diffusivity is present rather general results due 
to Avellaneda, Majda and Vergassola [23] show that, if 
the infrared contributions to v(x, t) are not too strong, 
standard diffusion occurs. 

Let us now discuss the well known results for the diffu- 
sive behavior of the standard map (a complete overview 
can be found in [22]), 



x(t + 1) = x(t) + y(t + 1) mod 2?r, 
y(t + 1) = y{t) + K sm(x(t)) . 



(16) 
(17) 



The Taylor argument, when applied to the y{t) compo- 
nent (in this work we always refer to the diffusion prop- 
erties of y(t)) of the standard map gives 



((xi(t) - x,(0)) 2 ) = 2D t + 11 <ftidt 2 <«i(x(t))u i (x(t))) 

Jo Jo 

rt rt2 oo 

= 2D t + 2 dt 2 dhCufo-h), (12) D E^ = ^ K 2^ in 2 x ^ + J2K 2 (smx(t)smx(0)}. (18) 



where we have assumed that (i>j (x(t))) = and we 
have introduced the correlation of Lagrangian velocities 
Cij(t) = (vi(x(t))vj (x(0))) . At large times, if the correla- 
tion decays sufficiently fast, the integral in (12) converges 
to an asymptotic value 



t=i 

At large K, the map (16-17) exhibits widespread stochas- 
ticity and to a good approximation consecutive angles x 
are decorrelated, thus one can neglect the second term in 
(18) to obtain the quasi-linear (or random phase approx- 
imation, RPA) result [22]: 



f 

Jo 



dtC u (t) = (v 2 )T L , 



(13) 



D*(K » 1) « Dql - 



K 2 



(19) 
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The above estimate is very crude: indeed it provides a 
good estimation of the diffusion coefficient only at very 
high K. Higher order corrections to the RPA approxima- 
tion can be obtained by means of the Fourier technique 
[22]. At order R- 1 ' 2 one obtains 



K 2 



8 ,v 57r > 



(20) 



This approximation is rather good apart for small K and 
around particular values of K. For K < K c w 0.972 
because of the presence of separating KAM tori there is 
not diffusion at all, Dy(K) = 0. On the other hand, 
at specific values of K (e.g. K w 6.9115) corresponding 
to the existence of ballistic solutions in the y direction, 
instead of the standard diffusion one observes an anoma- 
lous transport with v > 1/2 [24]. 



interaction with particles, evolve according to the gener- 
ating function 



S f = 9J' + J Lu{J')dJ' . 
according to which 



J' 



+ w(J'), 

J. 



(24) 



(25) 



That is, in the absence of coupling the amplitude and 
frequency of the mean field are constant. 

The self-consistent coupling between the particles and 
the mean field is specified by two functions / and g in 
the interaction generating function 



N 



Si = g( J') £ /(*„ - 9) . 



(26) 



III. SELF-CONSISTENT MAP MODELS 

In this section we introduce the symplectic map models 
that we propose for studying diffusion in self-consistent 
systems. As discussed in Section I, these maps consists 
of ensembles of globally coupled standard maps. 

The definition of the maps and the coupling is guided 
by the well-known fact that if (q, p) denotes the canonical 
conjugate coordinates of a Hamiltonian system at time 
n, then the transformation (q, p) — ► (q',p') given by 



q = 



dS 
dp' 



dS_ 
dq 



(21) 



defines a symplectic map with generating function S = 
5(q,p') [22]. 

The generating functions of the models proposed here 
have the generic form: 



S — S p + S f + Si 



(22) 



where S p defines the uncoupled evolution of the particles, 
Sf defines the uncoupled evolution of the mean field, and 
Si defines the particles-mean field interaction. 

For S p we assume the standard map generating func- 
tion 



N 



n=l 



X nV'„ + Vn + K n COS X n 



(23) 



where r„, and K n are constants, the index n labels the 
particles, and we have used the notation x(t) — x and 
x(t + 1) = x' . In S p the particles are uncoupled and each 
one follows independently a standard map dynamics. 

In order to preserve the symplectic structure of the sys- 
tem, the field is represented by two conjugate variables 
a phase 9 and an amplitude J which, in the absence of 



Based on this generating function, we will consider two 
models, one introducing a coupling only through the 
phase of the mean field, and another introducing a cou- 
pling through the phase and the amplitude of the mean 
field. 



A. Phase coupling 

In this case, it is assumed that 

g = e , f = sm(x n -8) w = Ct J' , (27) 

and r„ = 1, K n = K for n = 1,2... TV; with e, Q, 
K constants. The complete equations of motion thus 
become 



x' n = x n + y' n mod 27r, 

y' n = y n + K sin x n + e cos(x n - 6) , 

e' = e + nj', 

J' = J - eJ2 n cos(x n - 9). 



(28) 



The parameter e measures the strength of the coupling. 
For e = we recover the situation discussed in the previ- 
ous Section for N independent particles. One may expect 
that, in general, the inclusion of a coupling between the 
different particles will change the diffusive behavior of 
the system. If there were no feedback from the particle 
variables x n to the field J (i.e. no self-consistency) one 
could argue that the effect of the field variable 9 in (28) 
would be the same as that of a noise and thus will destroy 
the correlations in the y n variables. In this case one ex- 
pects that the deviations of the diffusion coefficient with 
respect to the quasilincar prediction discussed in the pre- 
vious Section would be strongly suppressed. Of course, 
this kind of argument can not be completely justified in 
presence of the full coupling in (28). 
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B. Amplitude and phase coupling 



IV. DIFFUSION IN PHASE COUPLED MAPS 



In this case, it is assumed that 

g = 2y/j>, f = T n cos(x n -6), w = -fi, (29) 

where f2 is constant and T n are N independent parame- 
ters. The dynamics in this case is determined by 



x n x n + Vnl^m 

y' n = Vn + K sin x„ - 2 Tn^fJ 7 sin(x„ - 

e' = e-n-jj. r„ cos (x„ - e), 

J' = J + 2 \[T' J2n=l r « S ™(Xn - 0) , 



(30) 



where, to simplify matters, we have assumed that the 
external field is such that K n /T n — K being K a con- 
stant. This globally coupled map was originally proposed 
in Ref. [10] as a symplectic discretization of the single 
wave model Hamiltonian system in Eq. (10). Compared 
with (28), in the map (30) J and 9 are both coupled to 
the particles and this leads to a self-consistent modifica- 
tion of the phase and the amplitude of the mean field. 
The map for J is implicit. However, rescaling the vari- 



ables y n 



T„ K, and defining 



N 



k = 2\fj , 7„ = 2 T„ , r] = 7" sin ( x n - 0) , 

71=1 

(31) 

the map can be written in a fully explicit form as [10] 



x n 
V'n 
K' 



x n + y' n , 

y n + K sin x n - k' sin (x n - 9) 
■? + r] 2 +r] , 
o _i_ J_ 9n 

" ~ «r Fir ■ 



(32) 



Let us study the diffusion properties in the phase cou- 
pled map (28). In particular, we will see that the effect 
of the mean field coupling is the randomization of the 
phase. That is, when coupling occurs only through the 
phase, self-consistency increases the stochasticity of the 
map, and the diffusive properties of the coupled maps 
are practically indistinguishable from the dynamics of a 
phase-randomized uncoupled map. 

In the limit e = the diffusion coefficient (for the y 
component), as a function of the parameter K, displays 
the complex behavior discussed in Section II. The first 
natural question is whether this complexity survives in 
the presence of coupling with a self-consistent mean field 
(i.e. for e > 0). 

In Fig. 1 we show Dy(K) (normalized with the RPA 
prediction K 2 /4) in the uncoupled case (s = 0) and for 
some other values of the coupling (e = 0.05 and e = 0.1). 
Numerically, Dy{K) is calculated by taking the large 
time limit of the expression: 



((y(t)-y(Q)) 2 ) 

2t 



(35) 



For comparisons, we also plot the value of Dy(K)/DQL 
with the higher order RPA corrections as given by (20). 
As one can see, for a generic large value of K the presence 
of a small coupling do not change the diffusion coefficient, 
apart for the values at which one can observe anomalous 
diffusion in the uncoupled limit (e = 0). Therefore, the 
first effect of the mean field is to remove the ballistic 
contributions to the dispersion. 



In this model, the mean field amplitude, k, plays the role 
of the standard map parameter (besides the K constant) 
which is self-consistently coupled to the particles through 
the order parameter ij. In the mean field particle dynam- 
ics one can define the total momentum of the system as 



Y 



JV 

^ ' "fn yn i 
n=l 



(33) 



where the first term on the right hand side represents 
the momentum of the mean field, and the second term 
the total momentum of the particles . This quantity is a 
constant of motion of (30), 



V = V . 



(34) 



This conservation law plays an important role in the dif- 
fusive properties of the system. 
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FIG. 1. Normalized diffusion coefficient vs. K for different 
values of the e parameter; squares for e = 0, circles e = 0.05 
and diamonds for e = 0.1. Solid line represents the theoretical 
prediction as given by the higher order RPA approximation. 
The number of particles is N = 60000, fl = 0.1 and the 
number of steps is 10000. 

Thus, let us try to focus on the case where the stan- 
dard map (s = 0) shows an anomalous behavior, for e.g. 
K = 6.9115, and see where the differences with the self- 
consistent map (e ^ 0) appear. The diffusion coefficient 
is an asymptotic quantity. For finite time, the evolution 
of Eq. (28) with small e maintains a memory of the be- 
havior at e = up to a time T(e) (saturation time). We 
dchne T(e) as the time at which the finite-time diffusion 
coefficient, as given by Eq. (35), is reasonably close to its 
asymptotic value. 

In Fig. 2 we show T(e), for a system with K = 6.9115, 
as a function of 1/ y/s. The approximately linear behavior 
in the plot indicates the dependence T(e) ~ exp(c/ v / e), 
with c an arbitrary constant. Let us note that the anoma- 
lous diffusion is mainly due to the presence of ballis- 
tic (non-chaotic) trajectories. Therefore the failing of 
anomalous transport can be seen as the recover of a 
generic statistical behavior. In this sense this is con- 
sistent with a scenario d la Nekhoroshev [22] . We also re- 
mark that, within the range of values investigated, there 
is not evidence of dependence on the size N as can be 
seen in Fig. 3 (again for K = 6.9115). It is worth men- 
tioning here that in Fig. 3 for e = one observes the 
typical anomalous behavior of the diffusivity, that is, 
D E oct 3 [24]. 
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FIG. 3. Dy /Dl against time for two different number of 
particles, K = 6.9115, and different values of e. Solid line is 
for N = 60000 and the dashed-line for N = 40000. 



The discussed results indicate that the main effect of 
the self-consistent field is to reduce the deviations from 
the statistical prediction. This is to be expected because 
in this mentioned before, the phase coupling 

leads to a randomization of the phase which is precisely 
what is assumed in the statistical arguments based on 
the random phase approximation. We check this state- 
ment by replacing the self-consistent field with an ex- 
ternal noise. We study a system of N particles whose 
evolution is now given by a time-dependent generating 
function 



N 



S({x' n }, {y„}, t) = Sv{x' n ,y n ) + e ]T f(x' n - v (t)) 



n=l 



n=l 



(36) 
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FIG. 2. Saturation time, T{e), vs. 1/y/e for K = 6.9115 
and N = 60000. 



where i] is a random process with the same statistical 
properties of the self-consistent 9, i.e., rj is a random 
number uniformly distributed in the interval [0, 2tt]. The 
result for Dy{K — 6.9115) /Dl against time is plotted in 
Fig. 4 for both the self-consistent and the random field. 
One observes that the random approximation is rather 
accurate. Therefore, taking into account the results for 
K large and the specific values of K where the standard 
map is anomalous, one can say that the dynamics for K 
large in the self-consistent map model (28) is equivalent 
to an effective standard map. 
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FIG. 4. D%(K = 6.9115)/D QL vs. time for different val- 
ues of e and N = 60000. The solid-line is for an external 
self-consistent field and the dashed-line is for an external ran- 
dom field. The random external field is generated by picking, 
at every time step, a random number regularly distributed in 
the interval [0, 2tt] . 

Most interesting, a priori, is the case of small values 
of K . As we have recalled in the previous Section, for 
K < 1 and e = there is no diffusion due to the pres- 
ence of KAM tori. However, as Fig. 5 shows, the phase 
coupled (e ^ 0) self-consistent map displays finite diffu- 
sion for arbitrary small values of K . This is, once again, 
a manifestation of self-consistent driven phase random- 
ization. 

The same scenario as for large values of K has been 
identified, that is, the external self-consistent field is 
equivalent, when diffusion properties are under study, to 
an external random field. This can be seen again in Fig- 
ure 5 where we also plot the diffusion coefficient obtained 
from the effective random map. At small K, one observes 
a weak dependence of D® on K, while Dy ~ e 2 /4, lead- 
ing to a finite diffusivity also for K = 0. As mentioned 
before, it is the breaking of the regular orbits of the stan- 
dard map (for K small) induced by the mean field, that 
allows the diffusion of particles. This is clearly seen in 
Figure (6) where it is shown some trajectories in the 
phase space for K = 1 for both the standard map and 
the coupled map model (28) with e = 0.1. 

In summary, for model (28) the coupling to an exter- 
nal self-consistent mean field is equivalent to the effect 
induced by an external random field. This result is valid 
also for K < 1 where the standard map shows barriers 
to transport. 
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FIG. 5. Dy vs K for different values of e. Here the num- 
ber of particles is N = 60000 and the final time is 10000. 
The solid-line is for an external self-consistent field and the 
dashed-line is for an external random field. The straight lines 
correspond to the lines Dy = e 2 /4 for e = 0.1 (dashed-dotted) 
and e = 0.05 (dashed-doubledotted) 




FIG. 6. Left, standard map phase space for K = 1, and, 
right, the trajectory of one particle for the system (28) with 
K = 1 and e = 0.1. 



V. DIFFUSION IN FULLY SELF-CONSISTENT 
MAPS 

In this section we study the diffusive properties of fully 
self-consistent maps. By this we mean maps coupled 
through the phase and the amplitude of the mean field. 
Our study will be based on the map (32) with = 0, 
and in the absence of external field (K = 0) since we are 
mostly interested in self-consistent effects. 

We consider an ensemble of particles with coordi- 
nates (xj, yj), j = 1, 2 . . . N, in the rectangle x £ (— w, n) 



7 



and y E (— 7r/2,7r/2). We will focus in the study of a 
Gaussian distributed active field 



7j = 7o exp 



< x j + yj) 

2a 2 



(37) 



For these nonuniform jj distributions it is useful to dis- 
tinguish between the particle variance cr^L, and the con- 
centration variance az,, defined as 



where 



([y - (v)p\ } P ■ 



i N 



(?)' 



(h/-(y)- 



i w 

TV E ^® 



(38) 



(39) 
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A. Subcritical diffusion 

In the first simulation we iterated the map with initial 
conditions k(1) — 0.8, and 9(1) = 0. The upper panel 
of Fig. 7 shows the initial 7 distribution. For t > 0, the 
scalar mixes and in the process modifies k and 9. In 
particular, in this case, as shown in Fig. 8, k oscillates 
in time around a mean value < K >— 0.966 (< . > is 
the temporal average) slightly below the critical value 
k c = 0.9716 for the destruction of KAM barriers and the 
onset of diffusion in the standard map. The phase (see 
also Fig. 8) decreases monotonically. 



0.03 y 




n^50 



FIG. 7. Phase space particle distribution in the fully 
self-consistent map (32) for a Gaussian distributed active 
scalar according to (37) with 2a 2 — 0.2 and 70 = 0.0269 and 
initial conditions k(1) = 0.8 and 0(1) = 0. The two panels 
show the particle distribution at the initial and final time, af- 
ter 137000 steps. The height of a vertical lines corresponds to 
the active scalar concentration jj of the j-th particle located 
at (x,y) = (xj, yj ). 



That is, the self-consistent coupling drives the system 
periodically between a diffusive regime with no KAM bar- 
riers (k > k c ) and a non-diffusive regime with KAM bar- 
riers (k < k c ). As shown in Fig. 9, this yields to diffusive 
particle transport in y, <j 2 y — 2Dt, even though on the 
average k is below the threshold for diffusion, that is 
there is subcritical diffusion. Because the peak of the 7 
distribution remains coherent, in this case there is no dif- 
fusion in the concentration, i.e. a 2 y = 0. As discussed in 
the previous Section, also in the phase coupled map (28) 
there is a subcritical diffusion regime. However, there 
is an important difference between these two cases be- 
cause in the phase coupled map the time evolution of 
the phase 9 exhibits random behavior while in the fully 
self-consistent map, as shown in Fig. 8, 9 has a regular 
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behavior. 

The oscillations of k are caused by the feedback effect 
of the active scalar trapped in the period-one island of the 
map (see Fig. 7). We describe the coherent part of the 
distribution, i.e. the part of the distribution trapped in 
the period-one island, as a macroparticle. The macropar- 
ticle representation is a sort of renormalization process 
in which a group of particles with different values of 7^ 
are replaced by one with an effective 7. The macropar- 
ticle concept provides a link between systems with large 
(or infinite) degrees of freedom and low dimensional sys- 
tems [16,12,25]. In the case considered here, at a given 
time t, the macroparticle rotates around the o-point of an 
effective standard map with coupling constant K = n(t) 
and phase 6(t). 

The oscillation period for a standard map of coupling 
parameter K can be estimated as [22] 



T = 



2tt 



arccos 



(1-f) 



(40) 




5 10 15 20 13718 13719 13720 

n n 



-5.337 





-5.338 

10 15 20 1.3718 

n " x 10" 

FIG. 8. Time evolution of k and 6 in the fully 
self-consistent map (32) for a Gaussian distributed active 
scalar according to Eq. (37) and initial conditions k(1) — 0.8 
and 6(1) — 0. The plots show the evolution in time 
windows at the beginning n € (0, 60) and at the end 
n £ (137000, 137060) of the run. The same periodic behavior 
is observed at intermediate times 
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FIG. 9. Subcritical diffusion in the fully self-consistent map 
(32) for a Gaussian distributed active scalar according to 
Eq. (37) and initial conditions k(1) = 0.8 and 6(1) = 0. Even 
though in this case, as shown in Fig. 8, the mean value of (k) 
is below the critical value for the destruction of all KAM bar- 
riers, the variance shows clear evidence of diffusive transport 
with 2Df = 0.0038. 

To use this result to calculate the rotation period of 
the self-consistent oscillation of k, note that according to 
the conservation of momentum in (33), 



N 



K 2 (t) = K 2 (l) + In Ml) - Vn{t)\ ■ 



(41) 



n=l 



In the macroparticle description, this relation can be 
written as 



k (t) — k (1) + Y mc Y mc 



(*)> 



(42) 



where T mc — J2 n 7« ^ s tne effective 7 e / / of the macropar- 
ticle and Y mc the y-coordinate of the macroparticle. Ac- 
cording to (42), the oscillation period of k equals the 
rotation period of the macroparticle which can be esti- 
mated using (40) with K = (k). For (k) w 1 this approx- 
imation gives a period T w 6 which is in good agreement 
with the numerical results (see Fig. 8). 



B. Self-consistent suppression of diffusion 

In the previous example, the constant rotation of the 
active scalar trapped in the period-one island gave rise to 
stationary oscillations of k and steady particle diffusion 
in y. However, this is not always the case, and it is pos- 
sible that diffusion is suppressed rather than maintained 
by self-consistent effects. As an example, consider the 
same initial conditions as before but with a smaller ini- 
tial value of the coupling parameter k, namely k(1) = 0.6. 
In this Fig. 10 shows, there is an initial regime 

in which k oscillates beyond k c and diffusive transport 
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is present with D — 0.0014. However, after a fraction of 
particles have migrated to regions of large y, k drops sys- 
tematically below k c and diffusion is suppressed. At this 
point the system enters in a transient subdiffusive regime 
leading to the eventual elimination of the diffusion. As 
shown in Fig. 10, and in more detail in Fig. 11, the sup- 
pression of the diffusion is accompanied by a damping 
of the coupling parameter n. Note that consistent with 
the estimation in (40), the period of oscillation remains 
constant T w 6. According to momentum conservation 
in Eq. (42), this damping can be viewed as a momentum 
transfer from the mean-field to the particles. This is rem- 
iniscent of the Landau damping mechanism in plasmas in 
which an energy transfer from the field to the particles 
leads to a colissionless damping of the field. 




x 10 



FIG. 10. Self-consistent suppression of diffusion in the map 
(32) for a Gaussian distributed active scalar according to 
Eq. (37) and initial conditions k(1) = 0.6 and 9(1) = 0. Up- 
per panel shows the time evolution of k and the lower panel 
the time evolution of the square of the particle variance. For 
n £ (l,6xl0 4 ), k reaches values above the threshold for KAM 
barriers destruction (k c ~ 1). For later times, the maximum 
of k drops systematically below n c and diffusion is suppressed. 
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FIG. 11. Damping of k during the self-consistent suppres- 
sion of diffusion shown in Fig. 10. The curve with larger 
amplitude shows the beginning of the time series (top panel 
of Fig. 10 with M — 25) and the curve with the smaller am- 
plitude shows the tail of the time series (top panel of Fig. 10 
with M = 14 x 10 4 ). 



C. Macroparticle instability and diffusion of 
concentration 

In the previous examples, the value of k was rela- 
tively small, and the macroparticle (i.e. the conglomer- 
ate of particles with the largest values of 7) remained 
coherent. This yields to either steady diffusion (see 
Fig. 9) or transient diffusion (see Fig. 10) of the particle 
distribution, but no diffusion of the concentration, i.e. 
o iy w 0. However, self-consistent effects can destabilize 
the macroparticle and the concentration field diffuses. To 
illustrate this, we consider the same Gaussian distribu- 
tion (Eq. (37) as before but with a larger initial value of n, 
namely k(1) — 3.3). In this Fig. 12 shows, up to 

n«5x 10 3 , k remains approximately constant, the phase 
9 n (panel (b)) decreases monotonically, the concentration 
variance er 7y does not grow and the particle variance a py 
exhibits standard diffusion. Around n w 5 x 10 3 there 
is a transition and k grows rapidly giving rise to a dif- 
fusion of the concentration and a jump in the particle 
diffusion. Figure 13 shows the active scalar distribution 
at two different times. 
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(a) (b) 




FIG. 12. Macroparticle instability and diffusion ol passive 
scalar concentration in the fully self-consistent map (32) for a 
Gaussian distributed active scalar according to Eq. (37) and 
initial conditions k(1) = 3.3 and 6(1) = 0. Panel (a) shows 
the time series of n n , (b) the phase 9 n , (c) the square of the 
concentration variance a^ y , and (d) the square of the particle 
variance a^ y . The lower and upper dashed lines have slopes 
equal to 1 and 0.62 respectively. Around n ~ 5 x 10 3 , the 
macroparticle looses coherence, K n exhibits large fluctuations, 
6 n drops to a constant value, and there is a jump in a^ y and 

2 

O p y. 




FIG. 13. Active scalar distribution during the macroparti- 
cle instability phase (panel (a)), and the diffusive phase (panel 
(b)) in the fully self-consistent map (32) for the run shown in 
Fig. 12. The initial condition corresponds to the Gaussian 
distributed profile in Eq. (37) shown in Fig. 7. In the plots, 
the height of a vertical line centered at (x,y) corresponds to 
the active scalar concentration jj of the j-th particle located 
at (x,y) = (xj,yj). 



D. Quasilinear diffusion 



In the phase and amplitude coupled map there is also 
a regime in which the phase is random and the map is 
equivalent to a random standard map. This leads to 
quasilinear diffusion as shown in the following. This 
regime is approached for initial k(1) such that, in the 
long-time limit, the time average value K e ff = l/(N t — 

Nq) J2n=N K ( n ) i s l ar g er than 1, being N t the final time, 
and No a proper time chosen to avoid the initial tran- 
sient behavior. Since the value of k depends on the it- 
eration time, the eddy-diffusivity is also time dependent, 
and we must define an effective diffusivity which is noth- 
ing but the time average of the instantaneous diffusion, 
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D eff = l/(N t - N ) J2n=N D y( n )- With theSe defini - 

tions, the quasi-linear approximation is recovered. This 
can be seen in Figure 14 by plotting the D e ff/DQL vs 
K e ff, and with solid line we plot the RPA approximation 
(20). 




FIG. 14. With circles the numerically obtained values for 
D e ff/DQL against K e ff. Solid line shows the RPA approxi- 
mation (20). 



active scalar loses coherence and this leads to a jump in 
the particle diffusion and to the diffusion of the concen- 
tration field; d) for large enough initial values of n there is 
widespread chaos and (with the exception of initial values 
close to accelerator models) self-consistency is shadowed 
by stochasticity leading to quasilincar difusion. The be- 
havior of the phase coupled map is in general different. 
In this case, we have shown that a) in the limit of large 
K, that is strong stochasticity of the flow, the external 
self-consistent field (at least for diffusion properties) is 
equivalent to a random field; b) the singular properties 
of the standard map, like the existence of ballistic modes 
giving rise to anomalous diffusion, are overcome by the 
external field; c) for K small, diffusion is different from 
zero, i.e. the external field breaks the barriers to trans- 
port present in the standard map. Moreover, again the 
external field is equivalent to a random driving field. 
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VI. CONCLUSIONS 



We have studied self-consistent diffusion in active sys- 
tem described by globally coupled symplectic maps. We 
focused our analysis in two systems: an ensemble of 
phase coupled maps, and an ensemble of maps cou- 
pled through the phase and the amplitude. The latter 
model is a symplectic discretization of the single wave 
model and as such represents a simplified description of 
self-consistent transport in plasmas and fluids. Numer- 
ical results obtained with this model indicate that self- 
consistency plays a critical role in the diffusive properties 
of the system. In particular, a) for small initial values 
of the standard map parameter, k(1) = 0.8, coherent 
oscillations of the active scalar give rise to periodic os- 
cilations of n above and below the threshold for barrier 
destruction (k c = 0.9716) leading to subcritical diffu- 
sion; b) for smaller initial values, k(1) = 0.6, there is a 
diffusive transient that eventually is suppressed by self- 
consistent effects; c) for larger values, k(1) = 3.3, the 
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